Turning light into a liquid via atomic coherence 
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We study a four level atomic system with electromagnetically induced transparency with giant 
X^^^ and x^^^ susceptibilities of opposite signs. This system would allow to obtain multidimensional 
solitons and light condensates with surface tension properties analogous to those of usual liquids. 



PACS numbers: 42.65.Tg, 42.50. Gy 

The applications of nonlinear optical media mostly rely 
on the adequate dependence of their refractive indices on 
the amplitude of light fields. It is well known that the 
figure of merit of a suitable material for practical devices 
includes a fast and strong response to the light field as 
well as low losses 1] which has motivated an active search 
for optical materials with the appropriate properties 0. 

On the other hand, a significant breakthrough in Quan- 
tum Optics has been the realization of giant optical non- 
linearities in ffases by means of atomic coherence and 
interference A technique that has attracted much 
attention is electromagnetically induced transparency 
(EIT) 0, H, 0, in which an opaque medium becomes 
transparent to a probe laser beam by the addition of an 
appropriate coupling laser beam. The adequate choice 
of an atomic level scheme and driving fields can yield to 
controlable nonlinearities with very interesting applica- 
tions in the design of nonlinear optical devices. This 
has been the basis for many studies on the resonant 
enhancement of nonlinear optical phenomena via EIT 
0, i, H, El [T5, la^. However, only a few of these 
works have investigated the formation of transverse soli- 
tons E2 mostly considering the role played by the 
giant Kerr nonlinearity. 

In this paper we study the optical properties of a sys- 
tem where atomic coherence can be used to control the 
dependence of the refractive index on the amplitude of 
the light field. Many novel nonlinear optical phenomena 
beyond the giant Kerr effect are described, the most in- 
teresting being the obtention of the so-called liquid light 
condensates 14], i. e. robust solitonic distributions of 
light with analogies to ordinary fiuid dropplets. 

We consider the propagation of a weak probe light 
field of frequency uOp in a medium composed of four-level 
atoms and a coupling light field of frequency uOc (see e.g. 
12]). A scheme of our system is shown in the inset of Fig. 
n In this kind of system, a coupling field of frequency uOc 
changes the level structure J. 5] and induces transparency 
for a probe beam of frequency uop. A second effect is the 
enhancement of the optical Kerr nonlinearity. 72, 73, 74 
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FIG. 1: Real (solid curve) and imaginary (dashed curve) 
parts of the susceptibility x given by Eq. for 7 = 30 MHz, 
72 — 10~^7,73 = 74 = O.OO67, Ai3 — A23 — 7, A24 — -1.57, 
= 7/2, as a function of \Q.p\^ = \iJL\^\Ep\^ /Ah^ . Inset: 
Schematic plot of the energy levels and optical couplings of 
the four- level atomic system. 



denote the decay rates of the atomic states and A 13, A23, 
A24 are light detunings. Direct electric-dipole transitions 
between the two ground states |1) and |2) are forbidden. 

Paraxial propagation along z of a probe laser beam 
through an optical medium, is given by: 



dEp 
dz 



En 



(1) 



kp = 27r/Ap and Ep are wave number and amplitude of 
the beam. The optical susceptibility x the rotating 
wave and adiabatic approximations for EIT in the pres- 
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FIG. 2: (Color online). Solutions of Eq. © for ^5 = 0.1 /im-\ 
1.2 /xm"^ and 1.6 fim~^ for our parameter set (values given 
in the text). 



ence of a coupling beam Ec takes the form [l^ : 



eon|r3pr4|A|: 



(2) 



where |^^p,cP = ImP |^p,cP/4^^ are squared Rabi fre- 
quencies and r2 = A23 - Ai3 - ^72, Ta = A13 + ^73, 
r4 = A24+Ai3-A23+i74 and A = B^^,\^\%\^ / {VlB), 
with 5 = r2 + |l^p|Vr4 + ll^clVrs, being 72, 73 and 74 
the decay rates of the atomic states and A13, A23 and 
A24 the light detunings. For our numerical examples to 
be presented later we choose an atomic density t] = 10^^ 
cm~^, an electric dipole moment ja = 3 x 10~^^ Cm (for 
alkalii atoms such as Rb or Ce, assuming for simplicity 

/^13 M23 /i24 = M), 72 = 10"^75 73 = 74 = O.OO67 

(with 7 = 30 MHz), A13 = A23 = 7, A24 = -1.57, 
= 7/2 and Xp = 800 nm. 

From Eq. we get the coefficients x^"^^ of the Taylor 
expansion of the susceptibility x = Xl^o X^^"^^^^ I^pI^"^ - 
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C = r2r3 + |i^cp, ^ = r3/r4-i. 

In Fig. n we show the real and imaginary parts of 
the susceptibility x as a function of the squared Rabi 




200 
z (lam) 



400 



FIG. 3: (Color online). Propagation of an eigenstate of Eq. 
(4) with /3 = 1.3 /im~^ in as medium with the full complex 
susceptibility given by Eq. (0). (a-d) 3D plots of \Ep(x,y, z)\^ 
for 2; = 0, 250, 350 and 450 /xm. (e) Beam width w{z). 



our parameter choice. As it can be appreciated from 
Re(x(^p)) (see Fig. QJ, the real part of the susceptibil- 
ity of the medium grows linearly with \Ep\'^ for low pow- 
ers (due to the effect of a positive Xr^) and decreases 
for high powers (due to a negative Xr^)i while the losses 
are comparatively small in this range. Thus, we have a 
balance of diffraction plus self-focusing for low field am- 
plitudes and self-defocusing for larger amplitudes. This 
type of competition is also found in media with the so- 
called nonlinearity of cubic-quintic type, i.e. those with 
a refractive index of the form n = no -\- n2\E\'^ -\- n4\E\^. 
These nonlinearities have attracted a lot of theoretical 
attention recently H El El EIS IHI because of their 
predicted ability, when 714 < 0, to prevent collapse of 
laser beams for sufficiently large powers, thus yielding to 
different stable two-dimensional light distributions 
The robustness of these light bullets has been recently 
connected with the formation of a liquid light condensate 
with surface tension properties similar to those of usual 
liquids These media are able to support stable vor- 
tex beams llil jiol [iQj and display interesting nonlinear 
phenomena [2l|. 

For our choice of parameters, using Eqs. 
and the relation n{Ep) 
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we obtain 

13^4 , 



-3.0154- lO-^'^m'^/V^, which 



frequency of the probe light \Q^p\^ = \ii\^\Ep\^ / Ah? ^ for ?/;£(r)e*^^e*^^, where (3 is the propagation constant. For 



are, respectively, ~ 10^^ and ~ 10^^ larger than those 
measured for usual nonlinear optical materials [l^ . 
These facts provide some analogies between our system 
and CQ media. However, the contribution of higher order 
and dissipative terms will be relevant for us. 

First we will construct stationary transverse self- 
trapped solutions of Eq. Q of the form: Ep{r^z) = 
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FIG. 4: (Color online). Evolution of two different eigenstates 
with = 1 for several propagation distances (in /im).For 
/3 = 1.5/im~^, the gain enlarges the radial size of the beam 
preserving the vortex (solid lines). Insets (a-c): Evolution of 
the eigenstate with (3 — 1.3/im~^ for 2: = 0, 180 and 190/im. 



^ 7^ the beam host a vortex of topological charge I. To 
this end we set Xi — ^ ^^^^ solve numerically the problem: 



1 d 



ar^ r ar ^ 



= 0, (4) 



with boundary conditions '^^(O) = and ^^^(oo) = 0. 

This gives us stationary beam shapes corresponding to 
different powers as a function of p. Let us first consider 
beams with i = 0. In Fig. [21 we show the results for 
P = 0.1 /im~^, 1.2 /im~^ and 1.6 /im~^. Low values of P 
yield to light distributions with quasi- Gaussian profiles. 
As p is incremented, the spatial shapes become narrower, 
but still keeping a Gaussian shape. For larger values of 
/3, the beam fiux grows rapidly and the peak intensity of 
the light distribution saturates due to the effect of a nega- 
tive n4, yielding to light distributions with almost super- 
Gaussian profiles. Due to the giant nonlinear response, 
these powers can be experimentally achieved by using 
mW continuous lasers provided the sources are highly 
stabilized in frequency (typically 1 MHz bandwidth). A 
warning is in order: for our parameter combination, the 
probe beam has a power smaller but close to that of the 
coupling beam thus a fully quantitative treatment should 
consider a vector extension of Eqs including the ef- 
fect of the probe beam on the coupling beam. Since this 
extension makes the analysis even more complex in this 
paper we restrict ourselves to the scalar model and the 
full vector model will be the subject of future research. 

We have used the eigenstates of the nondissipative case 
as input conditions for propagation in a medium with the 
full complex susceptibility of Eq.((2|) (i.e. including the 
imaginary part of x). We observe that for p = 1.2 /im~^ 
the eigenstate keeps its shape while propagating in such a 
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FIG. 5: (Color online). Coalescence of two eigenstates with 
P = 1.6 /im~^ launched in parallel with the same phases. 



medium. In this situation the amount of energy pumped 
into the soliton and taken out by the linear and nonlinear 
gain and dissipation respectively achieves an equilibrium 
[2^. Eigenstates with P < 1.2 /im~-^ tend to spread 
during propagation since they do not achieve the critical 
power for the formation of a soliton, while those with 
P > 1.2 /im~^ keep their peak amplitude and increase 
their width during propagation as it can be seen in Fig. 
I2I In this situation the energy available in the medium 
is the responsible for the broadening of the liquid light 
droplet in a similar way to the process of growth of a fiuid 
droplet in a supersaturated atmosphere. This means that 
although small, nonconservative effects play an impor- 
tant role in the propagation of wavepackets for this set 
of parameters. We have found numerically with a very 
high accuracy that the radius of the light droplet varies 
as R{z) ~ z^, which is faster than typical fluid droplet 
growth phenomena 23] or models similar to C-Q ones 
such as the Ginzburg-Landau equations Eij. In both 
cases the growth takes the form z^ with q < 1. Even 
pure diffractive propagation leads an exponent q = 1, 
which is smaller than the one observed in our system. 

Next, we have constructed eigenstates of Eq. with 
£ = 1. Eigenstates with P < 1.3 /im are unstable under 
propagation due to the presence of gain and the vortex 
breaks into two fundamental beams as shown in Fig. ^a- 
c). However, for P > 1.3 jam the beams reach a critical 
value of the energy SO that the liquid light condensate 
is formed and the surface tension is able to sustain the 
vortex within the beam. Thus, the effect of gain is to 
enlarge the radial size of the beam without destroying the 
vortex. As the beam propagates the width of the beam 
surrounding the vortex increases while keeping the peak 
density constant, which again resembles the growth of an 
incompressible fluid. It is remarkable that the maximum 
densities in Figs. El and ^ are very similar. 

To study the robustness of these solitons we have made 
a series of numerical experiments with collisions of dif- 
ferent beams (we show results for £ = 0). First we have 
launched initially parallel beams in phase corresponding 
to eigenstates with P = 1.6 /im~^. Their mutual effective 
interaction, as it happens with solitons of the ID nonlin- 
ear Schrodinger equation, is attractive and leads to their 
fusion and subsequent trasverse oscillations of the new 
bound state as it can be seen in Fig. [HI In a different 
series of numerical experiments we have studied the colli- 
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FIG. 6: (Color online). Collision of two hipergaussian beams 
launched with opposite phases, amplitudes A — 890 V/m and 
width w — \ /im, separated 4 /xm and being one of them 
slightly displaced in the y-axis. The initial angle is 0.25 rad. 
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FIG. 7: (Color online). Total reflection at a nonlinear-linear 
interface of an eigenstate with /3 = 1.6 /im~^ and an incidence 
angle of 0.017 rad, at z = 0, 22, 34 and 56 /xm. 



sions of the same beams launched initially with opposite 
phases as shown in Fig. El As it can be seen, both beams 
survive the collision behaving as light droplets, although 
slight transverse oscillations of the beams are observed 
after the collision by excitation of surface modes. Both 
phenomena arise in fluid droplets collisions [2^. 

Finally, we have launched an eigenstate with /3 = 
1.6 \imr^ with an incidence angle of 0.017 rad against 
the frontier between a nonlinear and a linear material. 
The results are shown in Fig. [3 where the breakup of 
the beam into smaller droplets is observed analogously 
to splashing fluid droplets |2^ . 

In conclusion, we have shown that the adequate choice 
of the parameters of a specific EIT scheme leads to a 
giant response for both and 77,4 (but different signs) 
and could allow to obtain stable two-dimensional liquid 
light condensates with surface tension properties similar 
to those of usual liquids. Our theoretical and computa- 
tional results could be the basis for real experiments in 
nonlinear optics with continuous mW lasers showing this 
phase transition and nice liquid-like properties of light. 
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